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0^ ■ Abstract. The stationary phase method is applied to diffusion by a potential barrier for 

an incoming wave packet with energies greater then the barrier height. It is observed 
that a direct application leads to paradoxical results. The correct solution, confirmed 
by numerical calculations is the creation of multiple peaks as a consequence of multiple 
l| reflections. Lessons concerning the use of the stationary phase method are drawn. 
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The stationary phase method (SPM), first introduced to physics by Stokes and Kelvin [1], provides 
an approximate way to calculate the maximum of an integral. It has in time become a standard tool 
in the armory, not only of physicists but biologists, economists etc. [2]. Below we shall briefly sketch 
the method. One of its main attractions is the apparent insignificance of details of the integrand with 
the exception of its phase. Already within its description, a series of limitations and assumptions arc 
made. While these are known to the experts they are often assumed implicitly and tested indirectly 
a posteriori by the success or otherwise of the results obtained. 

Recently much interest in the physics community has been stirred by the results of this method 
applied to tunnelling times in a potential barrier [3] . This has resulted in predictions of superluminal 
velocities, or more precisely to tunnelling times which, in the so-called opaque limit, are independent 
of the barrier length. Now, while not addressing this question directly in this paper, we investigate 
what we consider a simpler but related problem: The (non-relativistic) diffusion of an incoming single 
wave packet with energy spectrum totally above the barrier height. We first show that a direct 
application of the SPM analogous to the tunneling case (energy spectrum below the barrier height) 
also leads to surprising, not to say, paradoxical results. We have then performed numerical calculations 
which clearly display secondary reflected and transmitted peaks. This stimulates the assumption of 
multiple reflections which when combined with the SPM yields excellent agreement with our numerical 
calculations. The primary lesson that we draw is that the SPM without additional knowledge such 
as the number of wave packets existing is ambiguous and whence meaningless. For diffusion problems 
the conservation of probabilities can in principle be used to eliminate this ambiguity. 
Consider a complex integral over an unspecified range of the form 

1= I F(k)dk= [ \F(k)\exp[i9(k)]dk , (1) 



for which |.F(A;)| has a single maximum within the range of integration at k = k . If 9(k) varies 
sufficiently smoothly within the interval where |F(fc)| is appreciable, we can expand 9{k) about the 
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point k = k in a Taylor series 



where 



6{k) = 9 + (k- k ) 6' a + 0[(fc - k ) 2 

d9(k) 



6(k ) and 6'o 



dk 



If the modulus of F(k) is sufficiently sharply peaked we can neglect the second and higher order terms 
in the above series. This allow us to approximate the integral in Eq.Q by 

X « exp[ie ] J \F(k)\ exp[i (k - k ) 9' ] d{k - k ) , (2) 

However, if 9' is large the function of k which is to be integrated oscillates rapidly and, consequently, 
this integral will be practically null. A significant contribution occurs only when, for appropriate 
values of any parameters within 9(k), 

9'o = . (3) 

In this study we consider a modulated plane wave and are interested in the configuration space 
wave function in one dimension x, 



tp(x,t)= / \F{k)\exp[i\(k)]exp[i(kx- Et)]dk 



(4) 



with E = k 2 /2m for non-relativistic quantum mechanics. The F(k) may be a gaussian or similar 
modulation function. The total phase is 



9(k; x,t) = kx- —t + X(k) , 
2m 



(5) 



and the condition 9' = then yields the x-t dependence of the maximum or peak of \ip(x,t)\. For 
example, when A(fc) = 0, we obtain the group velocity result for a free wave packet 



x = — t . 

m 



(0) 



The existence of a X(k) produces a time or space shift 

k 



x = — t - A' 

m 



k„ ( _ m\' 

kg 



(t - At) 



It is exactly this type of analysis which leads to a delay time in the reflection of an incoming wave 
packet impacting upon a step potential when the momentum or energy spectrum is totally contained 
below the step height [5]. A similar analysis has been used for tunnelling times [6]. 

The standard procedure in these one dimensional potential problems is to find the stationary (but 
not normalizablc) plane wave solutions with the appropriate continuity conditions (see below) and 
then pass to a normalized wave packet by means of a modulating function. While the plane waves 
exist at all times in an infinite range of x, the wave packet is predicted by the SPM to exist for the 
incoming wave for say t < (with appropriate chosen time origin) while the reflected wave and other 
waves exist only for t > 0. Around t = we will have interference effects, due to the simultaneous 
presence of both incoming and reflected wave, and for the below barrier case we also have (over this 
transitory period) a wave function within the classically forbidden barrier region. 

In the following figure, we show the potential barrier 
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divided into three regions I (x < 0), II (0 < x < I) and III (x > I). The dotted line indicates the mean 
energy of the incoming wave $,„<;, 

^indx, t) = Jg(k) exp[i(fcx - Et)} dfc , (7) 

V2mVo 

with g(k) a truncated gaussian or similar, peaked at k {E a = k 2 a /2m). Truncation is needed, at least 
for small k as indicated in Eq.Q, since we wish to avoid any tunnelling phenomena. The x-dcpendcncc 
of the plane wave solutions in the three regions are given by 

Region I: x < , cxp[ifcx] + R(k) exp[— ikx] [k = \/2mE ] , 

Region II: 0<x<^, A(k) cxp[i<jx] + B(k) cxp[-iqx] [q = y / 2m(E - V ) ] , (8) 

Region III: I < x , T{k) exp[ifcx] . 

R(k) and T(k) are the reflected and transmitted amplitudes respectively. The coefficients A(k) and 
B{k) are the right and left going amplitudes in region II. All amplitudes are to be modulated by the 
function g(k) eventually. Continuity of ^f(x,t) and its derivative at x = and x = I determines the 
coefficients A, B, R and T, 

A{k) = k(k + q) cxp[iA(fc) - iql ] / V(k) , B(k) = k(q - k) exp[iA(fc) +iql}/ V{k) , , , 

R(k) = {k 2 -q 2 ) sm[ql] exp[i\(k) - i |] / V(k) , T(k) = 2kq exp[iA(fc) - ikl] / V(k) , ^' 

where 

V(k) = {Ak 2 q 2 + (k 2 - q 2 ) 2 sm 2 [ql}}^ and X(k) = arctan {(fc 2 + q 2 ) tw[ql } / 2kq} . 

To apply the SPM in what we would call the naive way, we must multiply each of the above amplitudes 
by the appropriate plane wave phases. For example, in the simplest case when g{k) is a real function 
we obtain 



(10) 



(11) 

The presence of the phase term X(k) implies a delay time in the reflected wave analogous to what 
happens for the step potential when E < V a . Since the phase of the incoming wave contains only the 
plane wave factors, i.e., it is devoid of a A(fc), the incoming peak reaches the barrier at x = at time 
t = (neglecting interference effects). For the reflected peak with the above expression for 9 R (k), we 
find the position of the peak of the reflected wave to be at 
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x = A'(fc ) - (k /m)t , 



(12) 



with 



A' (fco) 



2 (fc 2 + q 2 ) k 2 ql - (fc 2 - q 2 ) 2 sin[<?2 ] cos[<?Z ] 

q 



4k 2 q 2 + (fc 2 - q 2 Y sin 2 [ql] 



(13) 



Note that only x < is physical in this result since the reflected wave, by definition, lies in region I. 

The above expression for the position of the reflected peak simplifies around the "resonance" values 
for fc (q ) where 

sin[g n Z] = , i.e. q„l — nn , 

with n a non-negative integer. Assuming therefore, for simplicity, a sharp spectrum for g(k) peaked 
at one of these resonance values 



2 <?o 



(14) 
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This predicts a delay time for the reflected wave given by 

(k 2 +q 2 )ml 



At r R es 



X res (k ) 



2 k a ql 



Now consider the corresponding "delay times" for the A, B and T waves. In particular, 



At r 



— W 

1 I res 
K 



(fco)-g'(fco)/] 



At: 



2k q- 



(15) 



(16) 



This is the delay time at x = 0. We can calculate a delay time even for the 5- wave at x = since it 
exists in region II. However, this is later than the time of the £>-wavc at x = I since it is left-moving 
i.e. as a consequence of the q(l — x) factor in Ea. (|10fl . Suffice it to say that it arrives at x = at a later 
time than the departure times of either the R or A wave peaks. For completeness the transmitted T 
wave packet has its peak at the start of region III, x = I, at the time 



t' 



< i 



Now already one may be somewhat surprised to note that the appearance of the transmitted wave 
coincides with that of the reflected wave. However, the above results become paradoxical as soon as 
one realizes that for the time interval from t = to t = At r ^ s < Af^ s this solution is devoid of any 
maximum (R, A, B or T). During this time, at least, we are clearly in contradiction with probability 
conservation since the incoming wave peak has disappeared at time t = 0. By choosing the wave packet 
dimensions small enough we can say that there is an interval of time in which the naive SPM says there 
are no significant amplitudes anywhere in x. Note however that this is only an heuristic argument 
since a peaked configuration space packet runs counter to the above resonance approximation (peaked 
momentum distribution). There are also other incongruities in this naive application of the SPM. If 
one recalls the well known step case with E > V , single peak reflection occurs instantaneously (zero 
delay time). One might expect that our results tend to this case in the limit I — > oo. This is not the 
case. It is also possible in some off-resonance cases to find negative "delay times". In these latter cases 
the maximum of the reflected wave and incoming wave would exist contemporaneously. This situation 
also implies problems with probability conservation. 

Numerical calculations automatically conserve probabilities, at least to within the numerical errors. 
So to understand what is happening we performed such calculations and an example of these is shown 
in Fig. 1, where a complex gaussian modulation function 



g(k) 



8^? 



exp 



a 2 (k-k ) 2 



exp[— ikx ] 



has been used. It is to be noted that the choice of including a phase factor in g{k) simply shifts all 
times by a constant mx /k at resonance. These figures display the wave function in the proximity 
of the barrier for suitably chosen times. One clearly sees in these figures the appearance of multiple 
peaks due to the two reflection points at x = and x = I. This observation suggested the following 
analysis and imposed the subsequent interpretation. 

The R, A, B and T amplitudes may be rewritten as series expansions by considering multiple 
reflections and transmission in the potential discontinuity points, 



R — ^ ' R n — Ri B.2 

71=1 
OO 

A = ^ ] A n = Ai 



I |^ ) 



k — q 



B = B n — _£>i 
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exp [2 iql] 
exp [2 iql ] 
exp [2 iql] 



(17) 
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with 



k — q 2 k 2k(q~k) m 4 kq , , 

= fc - ^ ' = aT+^ ' = W cxp[2 ^ ] ' Tl = ¥TW cxp[l (9 " fc) ' ] ' 

R 2 = ^-A 1 B 1 , (18) 
A; 

Rn+2 ^4«+i B n j_ 1 T n +! 

-^ro+i -^n 



cxp[2 igZ ] 



These sums reproduce exactly the expressions in Eq.@. In this form the interpretation is easy. i? x 
represents the first reflected wave (it has no time delay since it is real). R 2 represents the second 
reflected wave. As a consequence of continuity, it is the sum, in region II, of the first left-going wave 
(Si) and the second right-going amplitude (A 2 ), i.e., 

R 2 = A 2 + B, = y A, B x . 

k 

This structure is that given by considering two "step functions" back-to-back. Thus at each interface 
the "reflected" and "transmitted" waves are instantaneous i.e. without any delay time. Indeed the 
SPM applied separately to each term in the above series expansion for R yields delay times which are 
integer multiples of 2 (dq/dE) l = 2 (m/q )l. This agrees perfectly with the fact that since the peak 
momentum in region II is q , the A and B waves have group velocities of q /m and hence transit times 
(one way) of (m/q )l. The first transmitted peak appears (according to this version of the SPM) after 
a time (m/q )l, in perfect accord with the above interpretation. 

Let us re-express what is happening. The incoming wave peak reaches the first potential discon- 
tinuity at x = 0. It instantaneously yields a first reflected peak (i? x ) and right-moving (A ± ) peak in 
region II. When this later wave packet reaches at time t = (m/q )l the second discontinuity at x = I, 
a part T x is transmitted into region III {x > I) while a part B 1 is turned back and eventually gives 
rise to the second reflected peak and so forth. Is this compatible with probability conservation? It is 
because of the following identity 

oo 

Y / (\Rn\ 2 + \r n n = i. (19) 

n=l 

This result is by no means obvious since it coexists with the well known result, from the plane wave 
analysis, 

oo oo 

\Rf + \T\> = \Y,Rn\ 2 + \Y, T "\ 2 = 1 ■ (20) 
n—l n—l 

In Fig. 2 we have re-plotted for various times the numerical calculations displayed in Fig. 1 and also the 
separate integral calculations based upon the above multiple pole model i.e. for particular R n (T n ). The 
latter wave packets arc represented by the curves. The former un-decomposed numerical calculations 
are plotted by various bullets. Agreement is excellent. 

In conclusion, the results of the SPM depend critically upon the manipulation of the amplitude 
prior to the application of the method. A posteriori this seems obvious. If we consider an amplitude 
say 

z(k;x,t) = |z|exp[ia] 
the SPM will yield one peak position for each given time. If we write the identity 

z — z 1 + z 2 , 

where z 1 = z — w and z 2 = w, and treat separately these terms, then the same approach will yield two 
peaks and so forth. The method is inherently ambiguous unless we know, by some other means, at 
least the number of separate peaks involved. Our above barrier analysis is simply a particular example 
of this ambiguity, for which we have presented a simple resolution, based upon multiple reflections, 
confirmed in detail by numerical calculations. 
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Fig. 1. The square of the amplitude modulus for different time frames. Only a fixed region in x close to the 
barrier is shown. Each figure should be multiplied by the adjacent factor, where it exists, to obtain the true 
curve. The parameters chosen for the plot are listed in the first frame. 



